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Lanchester-type combat attrition process. The attrition 
process is analysed in detail and modelled as a so-called 
renewal process in which times between casualties are 
considered to be independently and identically distributed 
random variables. Other random variables that can be 

considered in the renewal process are examined, and the 
distributions of these random variables are determined in 
order to study the behaviour of attrition process. Examples 
with specific distribution functions are given for better 
understanding. Computer simulation is generated and compared 
with the attrition model developed. The total casualty 
occurence by total force is also discussed, through pooling 
of the single renewal processes. The total casualty occu- 
rence is shown to be a Poisson process and times between 
casualties to be approximately exponentially distributed for 
large numbers of combatants. 
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I. INTRODUCTION 



Today, deter ainistic-differential-eguation models ( that 
are commonly called Lanchester-type combat models) are 
widely used in military, operations research for defense 
planning purposes, especially as concerns land combat. In 
combat analysis, such questions frequently arise as 'Who 
will win the battle ?, How long will the battle last ?, How 
many survivors will the winner have ?, and How do the force 
levels change over time ?',etc. There have been a lot of 
studies to estimate attrition in battle to answer such ques- 
tions as these. 

In 1914 ,F. N. Lanchester [Ref. 1] first proposed a mathe- 
matical model of combat in order to justify the principle of 
concentration of forces under "modern" conditions. He 
believed that in "ancient" times war was a series of 
one-to-one duels between men engaged in hand-to-hand combat 
using sword, axe, lance, etc. So tnere would be no advan- 
tage gained from concentration of forces. Lanchester 
reasoned that modern technology had changed the nature of 
this "ancient" warfare and now there was a decided advantage 
to "concentration of forces". 1 He hypothized that each 
side's casualty rate was proportional to the number of 
opposing firers. Implicit ( but not yet explicit) in 
Lanchester' s original work is the concept of an "attrition- 
rate coefficient" as the rate at which a single firer kills 
a particular enemy target type. Furthermore, Lanchester 



1 Helmtold notes that the analysis of the data from 92 
land combat battles suggests that victory in cattle is 
primarily determined by the factors other than numerical 
superiority and challenges the ability of any model of 
combat which concentrates almost exclusively on numerical 
force size to yield a practically useful predictor of 
victory in battle. 
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tnought that the number of survivors were the critical 
factor deciding the outcome of the battle in "modern" 
warfare. 

Since Lanchester's original efforts to describe the 
dynamics of combat mathematically, many analysts - in an 
attempt to add reality to the combat description - have 
extended the theory to include additional factors. Many of 
these extentions were described by Dolansky [Ref. 2]. In 
his paper, Dolansky pointed out that the Lanchester 
attrition-rate coefficients were aard to determine for 
particular weapons, and, accordingly, coefficients of known 
magnitude had been assumed in most models. Thus, Dolansky 
concluded that the use of Lanchester-type models for 
prediction of battle results had been hampered by this 
inability to predict numerical values for attrition-rate 
coefficients. However, even if we may be in doubt as to the 
proper form of attrition relations, we may at least be able 
to make pragmatic assertions such as, ’If the Lanchester 
attrition is used for predictions, then by measuring certain 
factors and performing certain operations a prediction of 
such-and-such an aspect of combat can be obtained. ’ So, 
although there has been a continuing discussion among mili- 
tary operations analysts about the merits of Lanchester-type 
attrition-rate, it is still generally accepted as the basic 
founding for all the combat models srudies and attrition- 
rate coefficient is still a core of these studies. 

3onder/5arf oot [Ref. 3] proposed that an attrition-rate 
coefficient be defined as the reciprocal of the expected 
time for an individual firer to kill an enemy target. A 
more thorough dicussion, however, of the justification for 
Bonder/Barefoot is to be found in Taylor [Ref. 4]. Barfoot 
[Ref. 5] suggested that the a ttrit ion p roces s be consi dere d 
as a renewal p roce ss . Thus, for a single Y firer, we will 
assume that each individual Y kills X targets according to 
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n-th kill, S ; 

n 

expected number cf kills in time t, ( t) , etc 



an attrition process in which times to kills are indepen- 
dently and identically distributed (i.i.d.) random variables 
(r.v). We will also assume that each hirer acts indepen- 
dently of any other hirer, and similarly for targets. Such 
an attrition process turns out to be a renewal process in 
the theory of stochastic processes. 

Once we introduce the theory of renewal process, we can 
consider many other random variables such as the time up to 

the number of kills in time t, N (t) ; the 

Thus this 

thesis starts with consideration of 3onde r/Taylor ' s concept 
of attrition-rate coefficient, and then studies the prob- 
ability distribution of the time between casualties 
according to the theory of renewal process. 

In this thesis, we will mainly focus on the fundamental 
Lanchester-type attrition paradigm which deals with homoge- 
neous forces, because understanding this basic paradigm is 
essential fcr extentions to further complex medals. 

This thesis is organized in the following fashion. 
First, we review a few relevant facts about the Lanchester 
attrition models. Then we reviev basic ideas of renewal 
process and develop the random variables whose study is of 
our concern, followed by the study of these random vari- 
ables. Then we will describe the combat attrition process by 
comparing the values of different random variables obtained 
through renewal process. Finally, the conclusion mainly 
contributes to the summary of the work and suggests further 
studies. 
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II. FUNDAMENTAL LAUCH E STER- T YPE ATTRITION PARA DIGM 



A. THE EASIC PARADIGM 



Let us consider combat between two homogeneous forces: a 
homogeneous X force opposed by a homogeneous Y force. We 
will focus on the f orce-on -force attrition process in tae 

homogeneous forces. The basic 



combat between these two 
Lanchester- type paradigm for "modern' 1 warfare 
the c asualt y r a te of suc h a hom ogeneo us f orce 



assumes that 
is directl y 



prop ort ional to the n umber of enemy forces, e . g. the X force 



te is 


given by 






dx 


dt 


- ay 


with x(0) 


= *0 ' 


orce 


casualty ra 


te is given 


by 


dy 


dt 


- bx 


witn y(0) 


II 

O 

> 



where ’a' denotes 

targets and is cal 

cient. Here x and 

0 

combatants (respec 

Many factors may 

the time being it 

functional depende 

problem of comput 

attrition-rate cce 

Two approaches 

to determine the v 

(1) A stati 

generated by 



the rate at which a single Y firer kills X 
led a Lan c neste r att rit io n- ra te coef f i- 

y denote the number of X and Y 
tively) at the beginning of the battle, 
affect the value of 'a 1 and 'b*, but for 
is not essential that we be explicit about 
nee of a and b. Let us next address the 
ing a reliable numerical values for such 
f f icients . 

have been developed in the United States 
alue of these coefficients as follows: 
stical estimate based on "combat" data 
a detailed Monte Carlo combat simulation. 
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(2) An analytical submodel of attrition process for the 
particular combination of firer and target type. 

In this thesis we will consider only the second approach. 
S. Bonder [Ref. 6] has called this approach the use of a 
free - standing or i n depende nt analytical model. The basic 
idea of this approach is to develop an analytical expression 
for each required kill rate by considering the single firer 
engaging a 'passive* target (i.e. one that doesn't fire 
back) and then tie all the attrition rates together in 
force-on-f orce combat with a Lanchester- type model. 

Barfoot [Sef. 5] in his study has suggested taki ng a 
L an ch ester att ritio n-ra te c oeffic ient as the reci pro cal of 
th e expected time for an i ndi vi dual firer to kill an enemy 
'tar get. * 2 TJithin context of the above homogeneous force 
Lanchester-type combat model, this means 

1 

a E[ T 1 , 

XY 



where T denotes a r.v. representing the time for a single 
Y firer to kill an X target and E[T] denotes expected time 
to kill. Taylor also supported Barfoot’ s idea in his recent 
study [Ref. 7] by comparing times between kills that were 
both exponentially distributed and non-exponentially 
distributed, and using previous equation to determine the 
expected time to Kill a target. 



2 However, the justir ication of this approach is not 
accepted and is apparently somewhat controversial. But in 
case of exponentially distributed times between kills, a 
deterministic Lanchester-type model may indeed be considered 
to yield the mean course of combat. Anyway, this approach is 
accepted as general guide for further study of attrition- 
rate coefficient. 
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B. DETERMINATION OF EXPECTED TIME TO KILL 



Bonder and Farrel £Ref. 8] have developed general meth- 
odology for determining the attrition-rate coefficient by 
determining the expected time to kill a target, S[T], for a 
wide spectrum of weapon system types. Morever, research 
since the mid-1960's has led to the development of several 
other methods for computing the expected time to kill a 
target. For present purposes, two methods will be of our 
interest, 

(1) Method based on sum of component event times, 
and 

(2) Method based on first-passage time in Semi-Markov 
process 

If we assume the special case of tactical interest: namely, 

the case of Mar kov - d epende n t fire, Bonder [Ref. 3] has shown 
that E[ T ] turns out to be as follows: 



(t +t_) (t +t J 

n r m r [ 1-P (H| H) ] 

Ef T ]=t + t -t + + {- 

a 1 h P(KJH) P (H| ) P (K | H) 



+ P(H|H) - ? } 

1 



where all symbols are defined in Table I This expession for 
E[ T ] holds for the following assumptions; 

(1) Markov-dependent fire with parameter ? , P ( H J H) , 
and P (HIM) 

(2) Geometric distrioution for the number of hits 
required for a kill with parameter P(K|H) 

In assessing the time factor to determine the 

attrition-rate coefficient, it turned out that target acqui- 
sition process and weapon system capaoility were most impor- 
tant factor from the formula. However, when we consider the 
time factor of the attrition process, target acquisition 
process would be the most important factor to study because 
it directly affects the required length of time to kill a 
target. So we will elaborate on target acquisition process 
for further use in attrition analysis. 
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TABLE I 

Variables for Expected Time to Kill 



time to acquire a target, t 

a 

time to fire first round after target acquired, t 
time to fire a round following a hit,t 

h 

time to fire a round following a miss,t 

m 

time of flight of the projectile, t 



1 



probability of a hit on first round, ? 



1 



probability of a hit on a round following a hit, P(H|H) 
probability of a hit on a round following a miss, P (H| M) 
probability of destroying a target given hit, P(KjH) 



C. TABGET ACQUISITION PROCESS 

Target acquisition is the initial step in engaging a 
target and a most deciding factor that affects the time to 
kill. An important distinction made in VEI’s Lanchester-type 
combat models is whether the target acquisition process of a 
single "typical" firer type is a seri al pr o cess or a 
process [Ref. 4]. The two modes for the target 
acquisition process considered by VRI's models ( including 
VECTOR-2 ) are as follows: 

(1) Serial Acquisition 

(2) Parallel Acquisition 
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Here a single firer using serial acquisition does not 
acquire targets while engaging another target/ say target A. 
When such a firer ceases to engage a target A ( due to kill, 
lost or any other reason), then he must acquire a new target 
all the way from the beginning. So it is assumed that he 
does not remember any acquisitions made prior to engaging 
the target A. On the other hand, a firer using parallel 
acquisition searchs for targets continuously, 'even while 
engaging target A and remembers those targets that have been 
acquired. Thus when such a firer finishes the engagement 
with target A, he can immediately snift engagement to the 
next target, provided chat such a target was acquired and 
processed during or before the engagement of target A. For 
both methods, the VRI models assume that a firer never 
engages a killed target again. 

"I • Serial A eg ui sit ion 

Now when we consider the basic Lanchester-t ype 
combat model in the previous section, a firer in serial 
process mode of target acquisition starts again from the 
beginning in searching a new target whenever the previous 
target has teen killed. Thus we assume that target acquisi- 
tion time is always needed to acquire a different target for 
engagement. Here it is assumed that the total-force kill 
rate is just the sin gle-f ir er-kill-ra te times the number of 
firers, e.g. 



dx 




since each firer is independent from any other firer. Also, 
in the case of a serial process for each firer, the expected 
time to kill a target is the same as has been shown in the 
equation for E[ T ] in section 3. A summary of attrition-rate 
coefficient results for serial acquisition is given in 
Table II . 



16 



TABLE II 

Attrition Process Osing Serial Acquisition Mode 



dx 

— = - ay 
dt 



a = 



E[T ] 
XY 



T = time for a Y firer to kill an X target 
XY 



(t + t ) (t +t ) 
h f D f [ 1-p (H| H) ] 

EfT] = t *t -t + + { +?(H|H)-P } 

a 1 h P (K| H) P (H| H) ? (K | H) 1 



2. Parallel Acquisition 



As we mentioned earlier, a firer in the parallel 
mode of acquisition continues to acquire new targets even 
while he is engaging a particular target. Cnee this target 
has been killed, then he can immediately shift fire to a new 
target which has been acquired already. In this case, 
total-force kill rate is given by the product of the kill 
rate of single firer against acquired targets and the 
expected number of firers who have already acquired one or 
more targets. 

The summary results for attrition process using 

parallel acquisition is given in table III . Note here that 

f is same as the orebabilitv that a single Y firer using 
xy 
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. 



parallel acquisition mode has available one or more acquired 

targets at which to fire at time t, because if he is firing 

we know that he already acquired the targets. Then it should 

also be noted that f v represents the expected number of Y 

xy 

firers who have already acquired one or more X targets. 



TABLE III 

Attrition Process Using Parallel Acquisition Mode 



dx 

— = - f a'y 
dt xy 



f : Prob[single Y finer using parallel acquisition 
xy is firing at X target at random time t } 



't 

f - 1 - exp { - x ' E (s) ds} 
xy - 0 x y 



1 

a* = :Y firer kill rate against acauired targe 

2[T* ] - 

XY 



I* = Time for a Y firer to kill an acquired X target 
XY (conditionai kill time) 



E : the rate at which a Y firer acquires X targets 
xy at time t when appropriate 



E[ T ’ ] 



(t +t ) 
h f 

t -t + 

1 h ? (K | H) 



(t +t ) 
m t 



P (HIM) 



i1 : P_(H|H)] 
1 P(K|H) 



♦P (H | H) -P } 
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using parallel a 
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delayed renewal 
to kill the fir 
acquisition but 
parallel acquisi 
already. But i 
firer who uses 
targets. 

Before c 
understand both 
and also Eonder' 
as the reciproca 
be recognized t 
primarily depend 



ay assume that it would take longer to kill 
firer usiag serial acquisition than a fxrer 
cquisition because of the time fraction of 
arget. This disfunction between these two 
s may be further studied as the case of a 
process. Here we assume that expected time 
st target would be obtained from the serial 
following time to kill would be from the 
tion, provided that tne targets are acquired 
n this thesis, we will only deal with the 
the same acquisition mode to attrit all 



ontinuing, it 
the basic La 
s approach of 
1 of expected 
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ent cn the acc 



is important that the reader 
nchester- type attrition model 
the attrition-rate coefficient 
time to kill a target. It must 
ted time to kill a target is 
uisition process. 
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III. ANALYSIS OF THE CO MBA T ATTRITION PRO CESS 
A. INTRODUCTION 

The lanchester attrition-rate coefficient is the rate at 
which a single firer kills a particular enemy target type in 
Lanchester-type combat. Development of technically-sound 
and scientifically-valid methodology for determining numer- 
ical values for Lanchester attrition-rate coefficient is an 
essential prereguisite for building a Lanchester-type combat 
model. The basic construct cf 3onder/3arfoot methodology was 
to take a Lanchester attrition-rate coefficient as the 
reciprocal cf the expected time for an individual firer to 
kill an enemy target as mentioned in the previous chapter. 
Taylor provided justification [Ref. 4] for taking Lanchester 
attrition-rate coefficient as the reciprocal of the expected 
time to kill. He started from the oasic hypothesis that 
combat is a complex random process, but it contains enough 
regularity that the appropriate Lanchester-type equations 
are a good approximation to the mean course of combat. It is 
clear that the casualty rate is equal to the reciprocal of 
the expected time for a force to inflict a casualty, when 
the times between casualties are exponentially distributed. 
However, in case the times between casualties are no longer 
exponentially distributed, Taylor used the suggestion made 
by Bonder and Barfoot. Bonder [Ref. 3] and Barfoot [Ref. 5] 
suggested defining the Lancnester attrition-rate coefficient 
as the expected rate which a single firer kills enemy 
targets. 

In the spirit of Bonder and 3arfoot, Taylor provided 
more rigorous justification [Ref. 7] for the Lanchester 
attrition-rate coefficient that does not assume an 
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exponential distribution for times between casualties. He 
considered the case where the initial force size of X and Y 
is large enough to insure a negligible probability that the 
battle is terminated before an attrition pattern could be 
established. He made no specific assumptions about the 
distribution of times between kills, but assumed that each 
individual Y force kills X targets according to an attrition 
process in which the times between kills are independently 
and identically distributed random variables. Thus in the 
parlance of the theory of stochastic processes, he said that 
such an attrition process is called a renewal process. 

Prior to this, attrition prediction has been difficult due 
to the inability to predict casualty patterns. 5y investi- 
gating the distribution of the times between casualties via 
renewal theory, combat attrition analysis becomes easier. 
Indeed, we know a great deal about the casualty pattern once 
we know the distribution of the interarrival time of kills. 
Is we consider more variables from the renewal process, we 
find more information about tne casualty pattern. 

Furthermore, we can track the development of the combat 
attrition more precisely by studying the probabilistic 
distribution function of those random variables that can be 
considered in the renewal process. 

We assume here that the reader is aware of the concept 
of a counting process. When we think about combat between X 
and Y forces, we can assume that time between kills by a 
single firer has some distribution F. Nov we are interested 
in observing the occurence of casualty and the number of 
casualties, N(t) , that have nappened in time interval (0,t]- 
Furthermore, when times between kills are i.i.d., we call 
this counting process as renewal process. 

In this chapter, the study of attrition analysis will 
include the following : 
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• The r.v., N (t) , total number of casualties by tine t , 
it’s distribution, and properties. 

• The expected number of casualties by time t. 

• Theorems discussing casualty occurence. 

• The limiting behavior of the casualty occurence. 

• The age and excessive life of the casualty occurence. 

• Examples with specific casualty distributions. 

• Application of renewal process for casualty estimation. 

• Total force attrition by pooling the single firer 
renewal processes. 

B. REVIEW OF RENEW AI PROCESS 

If the sequence of nonnegative random variables {X ,X , 
...,} are i.i.'d., then the counting process (N(t), t>0 } is 
said to be a renewal process (or often called an ordinary 
renewal process) . Thus, a renewal process is a counting 

process such that the time until the first event has seme 

distribution F, the time between the first and the second 
event has, independently of the time of the first event, 
the same distribution F, and so on. So when a event occurs, 
we say that a renewal has taken place. 

An example of a renewal process, let us suppose that we 
have an infinite supply of lightbulbs whose lifetimes are 
i.i.d.. Suppose also that we use a single lightbulb at a 
time and when it fails we immediately replace it with a new 
one. Under these conditions, {N(t) , t>0 } is a renewal 

process when N (t) represents the numoer of lightbulbs that 

have failed by time t. Accordingly, if we assume that this 
failure is a kill by a single firer - X or Y -, then killing 
process also can be considered as a renewal process. Thus, 
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renewal (or event or arrival) will be synonymous vitn kill 

or casualty occurence in the rest of this thesis. Tie now 

assume that casu al ty occur s nee is a k ind of renewal process . 

If, in particular, the distribution of interarrival time of 

kill is exponential with probability distribution function 
x 

e , then this renewal process turns out to be a Poisson 
proc ess with rate, . [Ref. 9] 

Thus we may assume that X casualty occurence due to a 
single I firer who uses the same acquisition mode to engage 
each X target can be explained by ordinary renewal process 
because we assume that times between casualties are all 
i.i.d. However, we often consider a counting process for 
which the first interarrival time has a different distribu- 
tion than the remaining ones. So, we may think of a Y firer 
who at first uses serial acquisition mode to engage an X 
target, tut once he has finished engaging an X target, then 
he may use parallel mode to engage the rest of the targets. 
So, in this case, we can assume that interarrival time of X 
casualty has different distribution between the first and 
the rest. 

formally, let [X , i > 1 } be a sequence of independent r.v. 

i 

with X having distribution G, and [X ,i>2} having distribu- 
1 i 

tion F which is different from G. Then counting process 

{N^(t),t>} is said to be a delayed renewal process. So, we 
can think of using delayed renewal process for different 
acquisition mode, but here in this thesis, we will only deal 
with the case of ordinary renewal process for simplicity. 

We now define for the renewal process a number of asso- 
ciated random variables whose study is the objective of this 
thesis. Table IV and Figure 3. 1 give the intuitive interpre- 
tation of the random variables obtained in renewal process. 

Here we can interpret each random variable as following: 

• X is the time between the (n-l)st and n-th X force kill 
n 
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TABLE IV 

Random Variables in Renewal Process 



i 



r. v 


interpretation 


X 

n 


time between the (n- 1) st and n-th renewal 
i.e. n-th interarrival time 


S 

n 


time up to n-th renewal 


N(t) 


total number ox renewals in (0, t] 


M(t) 


expected number of renewals in (0, t] 
known as renewal function 


2 (t) 


length of time measured from the last 
renewal until a given time t , called as 
age or backward recurrence time 


Y (t) 


length of time measured from time 
t until the next renewal, called as 
excess life or forward recurrence time 


S (t) 


length of renewal time at time t, 
called as spread 



by a single Y firer. 

• S is the time up to n-th kill. 

n 

• N (t) is the total number of kills in (0,tj. 
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1 




< S ( t) > 



— Z(t) — > 


< Y (t)--> 

time 




— — — — — — x 


1 


: S 


N(t) 


N (t) +1 



Figure 3.1 Renewal Process. 

• M (t) is average number ox kills in (0,t], 

• Z (t) is the time between the last kill and a given time t. 

• Y (t) is the length of time from time t up to the 

next kill. 

• S (t) is the length of lifetime of X force soldier at 

time t. 

In this thesis, we will only deal with X force casualty 
by Y firer to avoid the confusion in notation. 

C. DISTRIBUTION OF 1GTAL NU3BER OF CASUALTIES 

Suppose that for an ordinary renewal process a single Y 
firer uses the parallel acquisition mode to engage passive X 
targets from the beginning of combat, when we assume that 

targets are already acquired. Then, the distribution of 
N (t ) , the total number of X force casualties by a single Y 
firer, can be obtained at least in theory by first noting 
the relationship that the numb er of kills by time t is 
greater than or equal to n if and only if the n- th k ill 
occu rs before or at time t . 

Formally 
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( 3 . 1 ) 



N (t) >n <==> S <t 
n 

From 3.1, we obtain 



?{N(t)=n} = P{N(t)>n) - P {N (t) >n + 1} 

= P {S <t} - P {S <t) 
n n+ 1 

Now, since the random variable X # s,(i>1), are independent 

i n 

and have a common distribution F, it follows that S = 2L X 

n i=1 i 

is distributed as F , which is n-fold convolution of F with 

n 

itself [Ref. 10]. Here, F is the probability distribution 
of time between casualties. Therefore we obtain 

P£N(t)=n] = F (t) - F (t) with F (t) = 1 

n n+ 1 0 

and 



P { N (t) <n} = P{S >t} = 1 - F (t) 

n n 



(3.2) 



Therefore the probability distribution of N (t) , which is 

total number of casualties by single Y firer, can be 

obtained explicitly fcr ail n. 

The simplest special case of equation 3.2 is obtained by 

taking the renewal process to be a Poisson process. 

Then S has the special Erlang distribution with n stages 
n 

because each F is exponentially distributed. Also, it is 
known from the previous argument that N(t) has a Po isso n 
dist ribution with mean t. 

B. EXPECTED SOMBER OF CAS0ALTI2S 



1 



Renewal Function 



Expected number of casualties, E[N(t) ], which is 
defined to be the mean number of casualties in the interval 
(0 , t ] , is known as the mean value or r en ewal fun c tion , M(t). 
From equation 3.1, the expected number of casualties in time 
t can be obtained as following: - — 
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-1 (t) 



E[»(t) J 






= 7 X ?[N(t)=k} 
k= 1 

S') 

=ZL — P (S (t) =k} 
n- 1 k=n 

fco 

= Z- P [N (t ) >n} 

n= 1 

£2 

= 2 _ ?{S <t} 
n= 1 n 

= X F (t) 
n= 1 n 



(3. 3) 



Eut here it is generally very difficult to find M (t) 
directly from the distribution F. Thus, we will employ the 
Lapl a ce Sti elt zes Tr ansfor m (L.S.T. ) [Ref. 9] to find .more 
conveniently an expression for M(t). We generally need 
three steps to find M ( t) . 



Step 1: 

From the definition of L.S.T. , we know tnat L.S.T. 
of a (t) is 

* f -st 

a(s) = e da ( t) 

0 

and also 

* f' 1 -st 

F (s) =' e dF(t) 

0 



Ste p 2: 

Then from equation 3.3 
* *2 * 

M (S) =-- F (S) 

n= 1 n 

* n 

= Z ’ F ( s) ] for s>0 
n=1 

The above is a sum of geometric series ana can be written as 
follows: 



* 

a (s) 



_L ( !L_ 

l-F^S) 



(3. 4) 
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or equivalently 



* 

F (s) 



H(3) 



* 

1 +M (s) 



Step 3: 

Since the L.S.T. of a function uniquely determines 
the function, the renewal function determines the distribu- 
tion of the icill time and therefore determines the prob- 
ability law of the casualty occurence. Specifically there is 
a one-to-one correspondence between interarrival time 
distribution F, and renewal function, H(t). So if we have 
the distribution F, then we can find M(s) from the equation 
3.4. But to find a(t),we have to transform K (s) to K(t) by 
inverse application of L.S.T. 

Let us try the example. Remember that this example 
will be applied all the way through this chapter. 



Example 3. 1 

Suppose that each interarrival time to kill X ,(i>1) 

i 

exponential distribution with . Tnen L.S.T. of F is 

* -st 

F (s) = e dF (t) 



has 



(** -st 
e 

- 0 



->t 

e dt 







(P + s) t 

dt 



7v 

A +s 






0 



(A+s) e 




A 

A+s 



if s>-'A 



Then from equation 3.4 
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. 




Vj^s) 

1-(-V(A+s)) 







s 



where 1/s was an L.S.T. of F{t) , which had value of t. 3 So 
M(t) = /\t + c (some constant) , but we know that -1(0) = 0, 

because there would be no casualty at time 0. So, finally 

M (t ) = At. 



that interarrival times of casualties are i. i.d. ana expo- 
nentially distributed with mean 30 minutes. Then 

a) Wnat is the distribution of total X casualties at 5 
hours later ? 

b) What is the probability that 3rd X kill was occured at 
time t = 4 ? 

c) What is the probability that at time t = 2, no X casu- 
alty occured ? 

To answer these questions we know tnat M (t) = 2t from above 
example because = 2, then number of X casualties in time t 
has Poisson distribution with mean 2t. 



3 Let F (t) = t for t>0. Then L.S.T. of F(t) is 



Let us be more specific about this example. Suppose 



F (t) 




s 



1 -st 

— e I 



> 0 



1/s for s>0 
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Then 



- 10 n 

e (10) 

a) P{N (5) =n } = - — 

n ! 



3 -8 

8 e 

b) P {N (4) =3} = . 0236 and 

3 ! 



0 -4 

4 e 

c) P {N (2) =0 } = ---- = .0183 



2 . Ben ewa 1 Bqua tion 

Now an integral equation for M(t) may be obtained by 
conditioning on the time of the first casualty occurence. 
Opon doing so, we obtain 

roo 

a (t) = . s; N (t) | X ^ = x] d F (x) (3.5) 

however 



E[ N(t) |X = x ] = ) 

1 v. 1 ♦ M (t-X) 



if x>t 
if x<t 



(3.6) 



for if the first casualty occurs at time x, x<t, then from 
this point cn the process starts over again, and thus the 
expected number of casualties in (o,t] is just 1 plus the 
expected number of casualties in time t-x from the beginning 
of an equivalent renewal process. Putting equation 3.6 in 
3.5 yields 

M (t) = \ |* : 1 ♦ M (t-x) ]dF (x ) 



= F(t) + 



M (t-x) dF (x) 



(3.7) 
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where equation 3.7 is known as the rene wal et ua t ion and F is 
a known function and M is an unknown function to ne deter- 
mined as a solution to the integral equation 3.7. Ross 
provided the solution for equation 3.7 [Ref. 11] as 
following and the solution is accepted true for every case. 

Ex pecte d n umbe r of casualties up to time t is 

M(t) = F(t) * f t F (t-x) da (x) (3.8) 

J v 

oo 

where H (x) = .El F (x) 
n= 1 n 

So it allows to find H(t) using only the known distribution 
function F, which is the distribution of time between casualt- 
ies , thus avoiding finding individual F . 

n 

3. Renewa l Dens it y 

Now let's define the renewal dens ity as m(t) = 
dd(t)/dt. Remember that H(t) is the expected number of casu- 
alties up to time t. So n(t) may be regarded as the expected 
number of casualties per unit time; more precisely, the 
average numiec of casualties during the time interval from t 
to t +h is 

( t + h 

M (t+h) - M(t) =\ m ( x) d x (3.9) 

Differentiating the renewal equation, we obtain the renewal 
equation for the density m { t) 

( t 

m (t) = f (t) f (t-x) m (x) dx for t>0 (3.10) 

So, this knowledge of the density f of interarrival time of 
casualties is sufficient to determine the density m(x) of 
casualty occurence. 
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E. LIMIT BEHAVIOUR CF CASUALTY OCCUREHCE 

Ive know that total number of casualties, N (t) , go to 
infinite as time goes to infinite, when we assume that force 
size is infinite. Eut it would be useful co know at what 

rate N(t) approaches infinity. That is, we would like to be 
able to say something about limit of N (t) /t as time goes to 
infinite. 

From the definition of r.v. we know that 

S < t < S 
N (t) N (t) + 1 

wnere S must be the time of the casualty after time t. 

N (t) + 1 

By refer to [Ref. 10 ] we know that 
N (t) 1 

— > as t --> Co (3. 11 ) 

t u 

where u = S[ X j = \ [ 1 - F (x) ]dx , wnich is the 

i > 0 

expected time between casualties. 

Here the function 1/u is often called as the rate of the 
rene w al pro ces s. Thus the average number of casualties per 
unit time converges to 1/u . Then, how about the expected 
average number of casualties per unit time? Is it true that 
M(t)/t also converges to 1/u as t — > ? This result, known 
as the e lement ary renewal theorem, will be stated without 
proof [ Ref. 11 ]. 

E lement ar y R enewal Theorem 

M (t) 1 

> as t --> £0 (3.12) 

t u 

Thus whatever the distribution F is, the average number of 
casu a lties up to time t is approximately (for large t) , 
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a (t) 



(3. 13) 



t 



a 



Also, the a vera ge number of casualties in the interval 
<t, t+h) f or h>0 is approximated for large t as following. 



M (t+h) - M(t) as t — > fro 

u 

for any distribution F with mean u. 



( 3 . 14 ) 



Here equation (3.14) is known as B lack w ell 's the o re m when F 
is not lattice. 

From the example 3.1, where casualty occurence time is 
exponentially distributed with mean 30 minutes, the average 
number of X casualties up to time t is approximated by t/u, 
which is 2t. Then after 2 hours of combat, we can assume 
that the average number of X casualties by single Y firer 
will be 2t, which is 4 casualties. 

Ano the r li miti nq r esul t which may prove to be useful in 
renewal process also concerns the r. v. N (t) . It is shown 
[Ref. 12] that N (t ) has an asymp toticall y n orm al di st ribu - 
tion w it h m ean t/u and var iance t A/u 3 . Ihus 



lim 

t~> Cv^ 



F{N (t) <y} 




_Y t/u 

~2 3 

\ t/u 



\ 



(3. 15) 



Suppose of the same example 3. 1 that casualty occurence mean 
time is 30 minutes. Then u = 0.5 = 2[ T ] = 1/^ . So A = 2. 
Also Var£T] = l/A = 1/4 = 0.25 hr' , which is 900 mini here. 
Let’s assume that we had a 10 hours of combat action. Then 
expected number of X casualties is E[H(10) ] = t/u = 10/0.5 = 
20 and var[N(10) ] = t = 10x0.25/(0.5 f= 20. 

Suppose of the question that how many X soldiers will be 
required in order to with probability of 95% that at least 
more than one X soldiers will survive after 10 hours of 
combat ? This question may be answered as following. 
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I 



0.95 < P £N ( 1 0 ) <x} 




x - 20 ^ 

Y'T0~~~ i 



-j- > x - 20_\ 
4.4 72 



Now from the normal table, 0.95 percentile of normal distri- 
bution is 1.645. So 1.645< (x - 20) /4. 472, where the answer 
turns out to be x>27.36 . So we can say that X force needs 
at least 28 soldiers as an initial force size to be 95* sure 
that at least more than one X soldier can survive after 10 
hours cf combat. 

There also exist a a lte rnative normal approxima ti on , 
when we assume that for integer n, 

?{N(t)<n} = P {S >t} where S = X 

n n i = 1 i 



from equation 3. 1. 

Then P{S >t} = 1 - P {S <t} 
n n 



T I t - nu , 

<r 'pr 



(3. 16) 



So now we assume that S has normal distribution with mean 

, n 

nu and variance nr . If we apply equation 3.16 to previous 
example , 



— / 10 - 0. 5x \ 

0. 95<P £H (10) <x} = 1 - 1 

• o.5 pr ■ 



— ,10 - 0. 5 x , 

Then ^ — < 0.0 5 

~ o. 5 nr ] 



10 - 0 . 5x 
or m-- < 

0.5 ^ x 



1. 645 



If we solve for x, it turns out to be x>28.83. So answer is 
almost same when we assume some rounding errors. 



F. AGE AND EXCESS LIFE DISTRIBUTION OF CASUALTY OCCURENCE 

As we mentioned earlier in this chapter, we can also 
consider such random variables as Y ( t) , which is the time 
from certain time t until next casualty occurence. 
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which means residual lifetime of a soldier in coahat from 
certain time t if he is still alive, and Z (t) , which is the 

time from t since the last casualty occurence (see Figure 
3 . 2 ). 



s 

N (t) 




s 

N (t) + 1 


— “X” 


i 




Z (t) 


Y(t) 






S(t) 





Figure 3.2 Age and Excess Life of Casualty Occurence. 



This means 



Y (t) = S - t 

N (t ) +1 



Z (t) = t - S 



N (t) 



where Y (t) is called as the exces s life or forward recur- 
rence time and Z (t) is called as a ge or b ack ward recur rence 
time. Also 

S (t) = Z (t) ♦ Y (t) 

= s -t+t-s 

N (t) + 1 N (t) 

= s - s 

N (t) + 1 N (t) 



N g) + 1 



I - 

=1 i i=1 i 



= X 



S (t) + 1 



(3. 17) 
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Here S(t) is often called as spread and S(t) is also r.v. 

because X is r.v. too. 

N (t) + 1 



Suppose of the Poisson process that X. is exponential, 
then E[ S (t) ] = E [ X { t ) ] + E[ Z (t) ] = 1/ A ♦ 1/\ X = 2/ A = 2S[X.]. 
So, S (t) can be assumed as the length of the lifetime of the 
X force soldier which is in combat at time t. 
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IV. THE SgP EfiPOSIT IOa OF CASUALTY OCCURENCES 
A. INTRODUCTION 

In the previous chapter, we have considered in detail 
the theory of casualty occurence through the ordinary 
renewal process of a single Y force soldier. We deal in this 
chapter more briefly with the superposi ti on of se vera l 
rene wal pro ces ses, which is tota l X f or ce c asualt y occ u rence 
behavio ur. When we assume multiple Y force soldiers who 

shoot at X force soldiers independently of each other, we 
may assume that all these Y soldiers kill X soldiers by 
independent casualty processes. So, we can say that there 
occur multiple independent ordinary renewal processes, 
simultaneous in time. 
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Figure 4. 1 Superposition of Casualty Occurences. 
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Suppose that we have ' p' independent casualty processes 

in operation simultaneously, all with the same P.p.F. of 

casualty occurence time. Consider the sequence of casualty 

occurences formed by pooling the individual processes. Fig 

4. 1 illustrates the special case, when p = 3. Formally, we 

can say that if {N { t) , i>0}, for i = 1 ,2, ... , k, are 

inde- 

i k pendent renewal 

processes and N (t) = II (t), then (N (t),t>0} 

p i=1 x p is 

called pooled process which is total X force casualty occu- 
rence behaviour by k Y force 

soldiers. 



B. SOME GENERAL PROPERTIES 



1 . Norma l A ppro ximatio n of Po ole d Output 



From casualty occurence processes in the previous 
chapter, we assumed that there exist a mean casualty occure- 

2 

nee time. So let u and (J~ be the mean casualty occurence 

i i 

time and corresponding variance by the i-th Y force soldier 
and let u(p) and g-(p) be the mean casualty occurence time 
and corresponding variance of the total casualty occurence by 
k of Y force soldiers. Then the question is how are these two 



val ues related each other ? 

From previous stuav, we know tnat 1) N (t) is approxima- 

i 

2 3 

tely normal with mean t/u and variance, <y t/u . Therefore 

i ' i i 

we can say that N (t) is also approximately normal with mean 

p 

k k 2 3 

2F. t/ui and variance rr~ t/u . 2) Also assume now that 

i= 1 i=1 i i 

either the pooled process is a renewal process or at least 
the analogous result holds, then approximately N (t) is norm- 



al with mean t/u (p) and variance ry(p)t/ u (p). 
So, from 1) and 2) , we must have 



t k 1 

= t H 

u (p) i= 1 u . (4. 1) 
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then u (p) = 



k 1 -1 

(XI ) which is the harmonic mean or the u , 

i= 1 u i 

i 



2 2 
(p) k i 

and t = t H t 

3 i* 1 3 

u (p) u. 

l 



2 3 k Ci 

then j~(p) = u (p) ,2T 

1 = 1 3 

u 

i 



jt 1 -3 Jc Gi 

( i = 1 u 3 

i u 

i 



(4.2) 



Also if we assume that u = u and h~~ = for all i = 1,2, ..,k 

i J i 

which means all casualty occurences are independent and identi- 
cal, then mean casualty occurence time and corresponding vari- 
ance of the total casualty behaviour by k of Y force are as 
follows ; 

1 -1 u 

u (p) = ( - k ) = - 

u k 



2 1 -3 T <T J~ 2 

(T (?)=(* k ) — - ic =---=(- ) (4.3) 

u 3 2 k 

u k 



Also total number of expected casualties, E[ N (t) ], can 
he obtained here by pooling individual casualty occurences 
as following; 

k 

E[N (t) ] = Z. E£H.(t) ] = k E[H. (t) ] = k M (t) (4.4) 

p 1=1 l l 

where M (t) is expected number of casualty from 
single firer, which is t/u for large t. 
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let's look at an example that deals with pooling of 
individual renewal processes. 



Example 4. 1 



Suppose of Y force artillery shells from several 
sources are bombording the same X target. Each source hurls 
shells at the target at a rate of 40 shells per hour. Assume 
that the interarrival times of shells from each single 
source are uniformly distributed over an interval (0, a]. 

Then the guestion is 1) What is the probability that less 
than 825 shells will be hurled at the target from a single 
source during a 20 hour period and 2) If there are 5 sources 
from which such attacks are launcned, what is the prob- 
ability that more than 4125 ( which is 5 times 825 ) shells 

are hurled at the target in a 20 hour period ? 

To answer the first question, let A = 40 shells/hour 
= 2/3 shells/min. Then u = E[ X ]= 1 - 3/2 min, when we 

assume pcisson process. But interarrival times are uniformly 
distributed over an interval (0,a], So, u = E[X] = (0 + 
a)/2= a/2, where a = 3 minutes. Then Var[ X] = (a - 0)~/12 = 
3/4 min' . Now let N (t) = number of shells hurled in (0, t] 
from a single source. Since N ( t) is normal with mean t/u and 
variance ty/u 3 from equation 3.15, for t = 20 hours = 1200 
minutes, the probability turns out to be as following; 



1200 

325 

pr I 3/2 ' — / S 25 - 800 \ 

P (N { 1 20 0) <3 25} = c =— 

~ V 1200 - 3/4 / £ 16.33 

nw 

= £(1.53) = 0.9370 
and for second question. 

Let N.{t) = number of shells hurled in (0 , t ] by i-th source 
where i = 1,2, ... , 5 . and 



5 

N (t) = H N (t) = total number of shells hurled in (0, t] 

p i=1 i 
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by 5 sources. 

Then {N (t) , t>0} is the oooled process for which mean 
p 

ir.terarriva 1 time is, 

u 3/2 3 

u (p) = - = — minutes. 

k 5 10 



and corresponding variance is 
2 t- , 2 3 . 2 

<r <?> * < -v* > = 

Since N (t) is normal with mean t/u (p) and variance y (p) t/(^-' (p) 
for t= *1200 minutes. 



4 125 - 



1200 



-r- I 3/10 \ 7tr/ 125 

? {N ( 1200) <4 125} = c£ =9 

p — 1200 3/100 I 36.51 



3/1 

(3/10) 






= ?t>{3.42) = .9997 



So, ? {N ( 1 200) >4 12 5} = 0.0003 
? 



2 . Po i sso n Process of Pooled Output 

Now let’s assume that individual casualty processes are 

probabilistically identical witn common underlying C.D.F. of 

casualty occurence time F (x) . Let F (x) be the C.D.F of 

p 

casualty occurence time of pooled process. Then the question 

is how is F (x) related to F (x) ? This question can be 

p 

answered by using the idea of excess lire. Let Y (t) be the 

i 

excess life of an X force soldier ana Y (t) be excess life 

of total X force. Then we know that Y (t) = min[ Y jt) , Y (t^ 

.... , Y ( t) ] (see Fig 5. 1) . 

k 

Then 
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Figure 4.2 Super position of Excess Life. 



k 

? {Y (t)>y) = ?{Y. (t) >y} 

P i 



(4. 5) 



From previous study, we know that P.D.F of E.D. for individual 
renewal processes is (1/u) [ 1 - F (x) ] and then ve can say that 

1 

[ 1 - F (x) ] is P.D.F of E.D- for the pooled casualty 

u (p) p 

process. Thus, 



P{Y(t)>y} 



P{Y (t) >y} 
? 





F (x) ]dx then. 




F (x) ]dx 
P 



k 

F (x) jdx] from equation 4.5 



If we differentiate with respect to y 



4 2 



1 

[1 - F 

u(P) P 



(y) ] = 



i 

u 
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k-1 1 

F(x) ]dx] (- - )[ 1 - F (v) ] 
u 



Here we know that u(p) = k/u, so u ’ s are all cancelled cut. 
Then , 



C.D.F. of pooled output is 



1 ( k- 1 

F (y) = 1 - [ - [1 - F(x)]dx] [1 - F(y)] (4.6) 

p a y 



For example, think of the exponential distribution of 

- X x 

casualty occurence. Then F(x) = 1 - e , and u = 1 /\ . 

Then C.D.F of pooled casualty occurence is 

-Xx k-1 -\y 

F (y) = 1 - [ \ [e ]dx] e 

p y 



-AY k-1 - > y 

= 1 - (e ) e 



-Xky 

= 1- e for v>0 



which follows that pooled output is a Poisson process with 



rate k where k is the number of Y force combatants. 



Here, probably the most interesting properties of 
the pooled output refer to the ’local’ behaviour where k is 
large. Khintchine [ Ref . 13] has proposed that in the limit 
the numbers of renewals in non-overlapping intervals follow 
independent Poisson distributions, thus showing that in the 
lim it the poo l e d ou tp ut is a Poisson process . Also his 
proof does not require each interarrivai time to be identi- 
cally distributed. In attrition analysis in combat, it 
appears to be the same way, which means that if a large 
number of independent (and identical) casualty occurences 
are pooled, then the total force casualty occurence is 
approximately a Poisson process [Sef. 9] [Ref. 13]. Thus, 
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we can say that po o 1 6 d casualty occ u rence by total Y comba - 
tant s c an be an proxi mate! b v P oisson process if we assume a 
lar T e num ber of Y forces . 

Thus, finally, from the basic theorem of Poisson 

process, we know that if (N (t) ,t>0} is Poisson process with 

rate >0, then the interarrival times (X ,i>1} are i.i.d. 

i 

exponential candom variable with rate A . Thus, from 
previous results, t ota l X casualty occurences is a Poisso n 
pr oc ess with rate k/\ , w here k is the number of Y force 
combata nts, and int erarri val time s of X casualti es can be 
consid ered to be exp onenti a lly distributed w ith pa ramete r k 

A 
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V. COMPARISON OF RE NEWA L PROCESS WITH COMPUTER SIMULATION 



In this chapter we will look at the results of seme 
typical combat attrition process ( Y firers with different 
number of passive X targets) obtained through the theoret- 
ical renewal process model which is developed in this thesis 
and computer simulation using random number generator. 

As we mentioned earlier, total number of expected casu- 
alties by n Y firers turned out to be n M(t) for large t 
from the renewal theory, where we assumed that all casualty 
occurences were independent and identical. But it is clear 
that these expected casualties can not be bigger than the 
number of X targets. Now we will compare this result with 
simulated output to see whether we can find any interesting 
facts such as how many X targets will be most appropriate to 
apply renewal theory to combat attrition ? or what force 
ratio will be most applicable to renewal process ?, etc. 

In order to facilitate the simulations, the assumptions 
were made that n firers shoot independently at m targets, 
and that the times to kill the targets are uniformly 
distributed between [0,1]. When a target is killed, all 
firers shooting at it are assumed to select a new target at 
random from those still surviving. 

The procedure of the computer simulation is as follows 
(see Appendix A) ; 

1. We generate the initial kill time for each firer by 
using pseudo random number generator, where we assume 
that every firer's kill time is uniformly distributed 
between [0,1]. 

2. Each target is indexed by integer number. 

3. We need to decide which target does the firer snoot 
at ? The target selection rule is such that first we 
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generate a random uniform number between [0,1], then 
we take the integer value of (random number times the 
number of targets ) plus one, which designates the 
index of target that is assigned to firer 1. Same 
procedures are repeated until all the firers are 
assigned a target, with some targets possibly hcving 
more than one firer. Thus at first time step random 
number for selecting the target is generated n 
times, because we have to assign a target to each 
firer. These targets are killed eventually as time 
proceeds. 

Now, the target that will be killed first is the target 

which is assigned to the firer who has the minimum time to 

kill ( = Tmin) among the firers, where this time to kill was 
already generated. So we now know which target is killed. 

After the first target is killed, the number of targets 
killed is registered as one and the first time step ends. 
Now, in the second time step, tne new kill times between 
[0,1] are generated again for all firers who were shooting 
at the first target. The firers who were still in the 
process of killing a target at the end of the first time 
step continue to engage the same target in the second time 
step, where time to kill in the second step is old kill time 

( generated value at first time step ) less Tmin. 

Cnee the second target is killed, the second time step 
ends and the number of targets killed is now registered as 
two. Now we proceed to the third time step, where same 
target selection rule and firing rule are applied until all 
the targets are killed. 

After all the targets are killed, each time step length 
( which is the same as the time to kill the target in each 
time step ) is stored, and this is one replication (see 
Appendix E) . To obtain the expected value, 40 replications 
of the simulation were conducted and 6 data sets 
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(n =3, 2 = 6 ), (n = 3,m=10), (n= 3, ®=14) , ana (n=10, m = 6 ) , (n=10, 

m=10), (n=10, m= 14) , were examined. After 40 replications , 

the data were evaluated to find the expected number of kills 
through order statistics. This expected numner of kills was 
plotted against time. The results for n=3 and n=10 are 
shown in Fig 5.1. 

The theoretical results can be obtained through two 
methods, which are discrete approximation and simulation. 
In this thesis, we used discrete approximation based on the 
renewal function modelled at equation 3.7. It's formula is 
as follows; 



M (t) = i-1 1 + a (t-j) ] P . 

3 = 1 3 

where P represents the probability distribution 
of kill-time. 

Thus, expected number of kills obtained from the renewal 
process model using discrete approximation (see appendix C 
for program and data, where left column refers to the kill 
time and right column refers to the expected number of kills) 
is alsc plotted against time in Fig 5.1. 

It was found through both cases that the expected number 
of kills by the renewal process model may only be justified 
through the assumption of an infinite number of targets 
because as m gets bigger, the expected kill value approaches 
the theoretical result. Also it was discovered that at the 
beginning of combat, the attrition process almost followed 
the theoretical renewal process model. 
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EXPECTED NUMBER OF KILLS EXPECTED NUMBER OF KILLS 



THEORETICAL VS SIMULATED RESULTS 



T 

i 



3 RRERS 



12 3 4 

TIME 



10 FWERS 





Figure 5.1 Expected Kills From Theoretical And Simulated. 
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VI. CONCLUSIONS 



a. summary 



Since Lanchester proposed his model, many studies have 
tried to model the combat attrition process mathematically. 
Recently, the main focus has been on calculation of numer- 
ical values for Lanchester attrition-rate coefficient. This 
numerical determination stands at tne heart of casualty 
assessment in such models. There are three reasons for the 
importance of this attrition-rate coefficient methodology. 
First, Lanche ster- t ype models are used in various U.S.Army 
and D.O.D. planning activities more widely than ever before. 
Second, Lanchester attrition-rate coefficient is a basic 
element of any Lanchester-t ype combat model, and that 
attr itior-rate coefficient reflects tne effective applica- 
tion of firepower. Finally, significant new developments 
have occured in methodology for developing more tactically 
realistic Lanchester attrition-rate coefficient and these 
important results have not been accessible to a very wide 
audience. In particular, a new approach for developing more 
realistic and mathematical Lanchester attrition-rate coeffi- 
cient, that of computing combat attrition by using a renewal 
process, was proposed by Taylor recently. This thesis began 
with the assumption that the combat attrition process was a 
kind of renewal process, and that times between casualties 
were i.i.d. Thus, considering the problem from the stand- 
point cf renewal theory, many new random variables were 
included and various distribution functions studied in order 
to understand better this type of 
paradigm. 



4 S 



casualty occurence 



This thesis only dealt with the case of homogeneous 
forces which were using the same acquisition mode to attrit 
the targets, which were passive X forces. Thus we could use 
an ordinary renewal process to study the behaviour of casu- 
alty occurences. We emphasized the probabilistic function 
of each r.v., so once we found the distribution of interar- 
rival times of casualties, casualty analysis could be done 
by plugging distribution functions into the eguation devel- 
oped. Finding the distribution of interarrival times of 
casualties was the critical factor in continuing the casu- 
alty analysis. 

The final conclusion is that total casu al tv o ccurence in 
co mbat can be assumed to be a Poisso n pr oces s , and thus 
times between c asu a lties are e xpo nential ly dist ri buted if ye 
assume a large n um ber of combatants. 

We have dealt here only with a passive X force ( i.e. one 
that does not shoot back ) , perhaps an unrealistic assump- 
tion. But still this study permits better understanding of 
combat attrition, through the use of probability theory, 
than was previously available. 



B. SUGGESTIONS FOR FURTHER STUDY 



In this thesis, we 
targets. A core complex for 
combat, since enemy targ 
attrition. 

In addition, a delayed 
for the case where differen 
the first shooting and for 
reinforcement of the fo 
Casualty analysis could th 
alternating renewal proces 
reinforcement. 



have only dealt with passive 
mulation may be required for real 
ets shoot back and cause firer 

renewal process could be applied 
t acquisition modes are used for 
each subsequent shootings. A 
rces could also be included, 
en be done using the theory of 
ses, once we know the rate of 
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APPENDIX A 

FORTRAN PROGRAMMING FOR SIMULATING T3E KILL TIME (N = 3) 



SAL U1(3) .02 (3) £ U3,U4,TMIN,TIME (3) , INTIME (40,14) 
NTEGER JIG1, F INDEX, L ,M, N, KIGT, IGT (14), RTGT (14) 



100 

105 



110 

10 



130 



125 



70 



135 
20 0 

300 



R 

INTEG 

DOUBLE PRECISION DSEED 
ESEED = 123457. DO C 
DO 200 K = 1,40 
M = 14 
N = 3 
L = 1 C 

CALL 3 GUBS ( DSEED, N, U 1 ) C 
DO 100 I = 1.N 
TIME (I) = U1(I) 

CONTINUE C 
DO 105 J = 1 , M 
RTGI(J) = J 
CONTINUE C 

DSEED = DSEED + 23. DO C 
CALL 3GUBS (DSEED, N, U2) C 
DO 110 I = 1 , N 

TGT(I) = I NT (U2 (I ) * FLOAT (M) ) + 1 

CONTINUE C 
CONTINUE 
TNI N = 2. 

DO 130 I = 1 , N 

I F (I MI N. L E . TIME (I ) ) 

T MIN = TIME (I) 

KIGT = TGT (I) 

FINDEX = I 
CCNIINUE C 

INTIME (K,L) = TMIN C 

DO 125 J = 1,M 



GO TO 130 



IF (RIGTjJ) .NE.KTGT) 



GO TO 125 



IT (M) 



JIGT 

RTGT(JTGT) = 

CONTINUE C 
M = M - 1 

I = L + 1C 
DO 135 I = 1 , N 

IF (KIGT. EQ. IGT (I) ) GO TO 70 
TIME (I) = TIME I - TMIN 
GO TO 13 5 

DSEED = DSEED + 1 2. DO 
CALL GGU3S (DSEED, 1,U3) 

TIME (I) = U3 

DSEED = DSEED + 3 . DO 

CALL GGU3S (DSEED, 1, J4) 

TGT(I) = RTGT(INT (U4*?LOAT(M) ) +1) 
CONTINUE C 
IF (M.3Z. 1) GC TO 10 
CONTINUE 

WRITE (6, 300) ( (INTIME ( X , L ) ,L=1, 14) ,K=1,40) 

FORMAT (IX, 14F5.2) 

END 
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SIMULATED DATA 



A? gENDIX 3 
FOB KILL TIME 



0 . 2607 
0 . 3914 
0 . 1156 
0. 4545 
0 . 4138 
0.0884 
0. 4117 
0. 0499 
0. 4637 
0 . 7428 
0 . 2326 
0 . 1105 
0. 0943 
0 . 5372 
0 . 0492 
0. 2496 
0 . 453 0 
0. 2916 
0 . 2299 
0.7527 
0. 2201 
0. 6057 
0.7823 
0 . 1133 
0. 003 5 
0 . 0174 
0. 6032 
0. 3186 
0 . 556 4 
0. 5706 
0 . 2102 
0. 019 1 
0.1570 
0 . 193 1 
0. 1372 
0. 0523 
0. 3024 
0 . 064 5 
0. 0220 
0 . 128 5 



0 . 5056 
0.3860 
0 . 2987 
0 . 0154 
0 . 1269 
0 . 0730 
0 . 3926 
0 . 24 10 
0 . 21 88 
C . 0736 
0 . 3066 
0 . 50 15 
0 . 1942 
0 . 0 1 04 
0 . 0 5 54 
0 . 0 2 27 
0 . 0027 
0 . 2342 
0 . 0 756 
0 . 0343 
0 . 3 036 
0 . 0538 
0 . 0 719 
0 . 6589 
0 . 0602 
0 . 56 69 
0 . 00 C 3 
3 . 04 92 
0 . 01 05 
0 . 3 844 
0 . 4294 
0 . 2929 
0 . 0 6 0 1 
0 . 4658 
0.41 83 
0 . 2950 
0 . 5260 
0 . 1 504 
0 . 4350 
0 . 1379 



0.2000 
0.1659 
0.01 87 
0.1312 
0.0678 
0.2945 
0.03 63 
0.3564 
0 . 1503 
0.36 29 
0.2037 
0.2316 
0.2346 
0.36 58 
0.1336 
0.1102 
0.1 524 
0.0265 
0.50 40 
0. 1628 
0 . 13 23 
0.05 16 
0.1470 
0.0144 
0.44 10 
0.1626 
0.26 92 
0.39 18 
0.0674 
0.35 04 
0.1155 
0.4283 
0.3111 
0.0251 
0. 1055 
0.4356 
0.0669 
0.0088 
0.2320 
0. 1063 



0. 4392 
0.3315 
0.391 1 

0. 2905 
0.1365 
0.0234 
0 . 0 o 97 
0.0033 
0 . 054 o 
0. 0375 
0.2982 
0 . 1160 
0.0342 
0 . 143 1 
0. 1750 
0.2093 
0 . 0436 
0.5775 
0. 1737 
0 . 0 063 
0. 5570 
0 . 333 1 
0.0156 
0.0712 
0 . 2291 
0. 2062 
0.0024 
0. 1603 
0.2037 
0 . 2022 
0 . 0393 
0 . 2152 
0 . 4439 
0 . 1568 
0 . 2 174 
0 . 0230 
0.0079 
0 . 0741 
0 . 0207 
0.0695 



(N = 3, 



0.1057 
0. 0303 
0. 0488 
0.1360 
0. 1 724 
0.7305 
0.0277 
0.1442 
0 . 5555 
0.0319 
0. 1050 
0.2916 
0 . 277 1 
0 . 0724 
0.5712 
0.2394 
0.0132 
0.3016 
0. 0637 
0. 3606 
0.0256 
0. 1 933 
0.0743 
0.1216 
0 . 1626 
0 . 1578 
0.0723 
0. 1857 
0.0134 
0.0417 
0.1705 
0 . 1427 
0. 0499 
0. 0375 
0.0335 
0. 1054 
0.1682 
0 . 553 1 
0.2634 
0 . 0050 



M = 6) 



0.0649 
0.1646 
0.4744 
0.3572 
0.0101 
0. 03 20 
0.0151 
0.2979 
0.0280 
0.0322 
0.1466 
0.4032 
0.0218 
0.1435 
0.0609 
0.1294 
0.6276 
0.1875 
0.5797 
0.0877 
0.3382 
0.0139 
0.0118 
0.2234 
0.2609 
0.0334 
0.2653 
0.0369 
0.1233 
0.1048 
0.4165 
0.4954 
0.0377 
0.2395 
0.2340 
0.0926 
0.0748 
0.0228 
0.2070 
0.5903 
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APPENDIX C 

DISCRETE APPROXIMATION FO RENEWAL FUNCTION 



$ JOB 

REAL PO, M 1 (401) ,11(41) , TIME (41) 

Ml (1) =0. 

PO = 1./100. 

DO 100 I = 1,400 

IF (I.GT. 1 00) GO TO 200 

Ml (1+1) = Ml (I) + (1. + Ml { I) ) * PO 

GO IO 100 

200 Ml (1+1) = Ml (I) ♦ (1. + Ml (I)) * PO 
* - (1. + M 1 (1-100) ) * PO 
100 CONTINUE 



DO 1 1 0 J = 1.41 
TIME (J) = 0.1* 


a, 


M { J) = Ml (1+ 10* 
WRITE (6,*) TIME 
110 CONTINUE 


STOP 




END 0 SENTRY 


0. 0000000 


0. 0000000 


0. 1000 00 0 


0. 1046219 


0. 2000000 


0. 2201895 


0. 300000 1 


0.3473478 


0. 4 00000 1 


0. 4833622 


0. 5000 00 1 


0.6446298 


0. 600000 1 


0. 3166940 


0. 7000002 


1. 0067590 


0. 3000002 


1. 21 67040 


0. 9000002 


1. 44861 60 


1 . oocoooo 


1. 7047900 


1. 0 9 999 90 


1. 8733930 


1. 1 999990 


2. 0587170 


1. 3 000000 


2. 2452650 


1. 3999990 


2. 4373700 


1. 5000000 


2. 6341 510 


1. 6 0 0000 0 


2. 83 44850 


1. 6999990 


3. 0369570 


1. 3000000 


3. 2393290 


1. 8999990 


3.4409610 


2. 0000000 


3. 6377740 


2. 1 00000 0 


3. 8320340 


2. 199999 0 


4. 0281400 


2. 3000000 


4.2256030 


2. 4000000 


4. 4239440 


2. 5 0 00000 


4. 6227140 


2. 6000000 


4. 8215160 


2. 6 9 99 9 9 0 


5. 0200490 


2. 8000000 


5. 2131320 


2. 9000000 


5. 41 57760 


3. 000000 0 


5. 6132290 


3. 1000000 


5. 81 09340 


3. 2000000 


6. 0089100 


3. 3000000 


6.2070170 


3. 4000000 


6. 4051470 


3. 5000000 


6. 6032300 


3. 6000000 


6. 80 12370 


3. 7000000 


6. 9991710 


3. 3000000 


7. 19 73 oOO 


3. 9000000 


7.3949500 


4. 0000000 


7. 5923790 
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